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ABSTRACT 

We  report  the  results  of  a  globally  coordinated  photometric  campaign  to  search  for  transits  by  the  P  ^  30  day  and 
P  ~  60  day  outer  planets  of  the  three-planet  system  orbiting  the  nearby  M  dwarf  G1 876.  These  two  planets  experience 
strong  mutual  perturbations,  which  necessitate  the  use  of  a  dynamical  (four-body)  model  to  compute  transit  ephe- 
merides  for  the  system.  Our  photometric  data  have  been  collected  from  published  archival  sources,  as  well  as  from  our 
photometric  campaigns  that  were  targeted  to  specific  transit  predictions.  Our  analysis  indicates  that  transits  by  planet  c 
(P  ^  30  days)  do  not  currently  occur,  in  concordance  with  the  best-fit  i  =  50°  coplanar  configuration  obtained  by 
dynamical  fits  to  the  most  recent  radial  velocity  data  for  the  system.  Transits  by  planet  b  (P  ~  60  day)  are  not  entirely 
ruled  out  by  our  observations,  but  our  data  indicate  that  it  is  very  unlikely  that  they  occur.  Our  experience  with  the 
G1  876  system  suggests  that  a  distributed  ground-based  network  of  small  telescopes  can  be  used  to  search  for 
transits  of  very  low  mass  M  stars  by  terrestrial-sized  planets. 

Subject  headings:  planetary  systems  —  planets  and  satellites:  general  —  stars:  individual  (Gliese  876) 


1.  INTRODUCTION 

The  G1  876  planetary  system  ranks  with  the  most  remarkable 
discoveries  that  have  emerged  from  the  first  decade  of  extrasolar 
planet  detections.  As  described  by  Rivera  et  al.  (2005),  Laughlin 
et  al.  (2005),  Butler  et  al.  (2001),  Marcy  et  al.  (1998,  2001),  and 
Delfosse  et  al.  (1998),  there  are  a  number  of  reasons  why  the 
G1 876  planets  are  extraordinary.  Planet  d  with  its  1 .974  day  orbit 
and  M  sin  (/)  =  5.9  Mq,  is  the  lowest  mass  exoplanet  known  to 
orbit  a  nearby  main-sequence  star  (Rivera  et  al.  2005).  Planets  c 
and  b  are  the  only  Jovian-mass  companions  known  to  orbit  an 
M  dwarf,  and  their  clear  participation  in  2 : 1  mean-motion  res¬ 
onance  gives  important  dynamical  clues  to  the  formation  and  evo¬ 
lution  of  the  system  (e.g.,  Lee  &  Peale  2002).  Furthermore,  the 
G1  876  planets  are  the  closest  exoplanets  that  have  been  reliably 
characterized.  Their  radial  velocity  amplitudes  induce  an  overall 
signal-to-noise  ratio  of  nearly  1 00,  and  they  have  been  observed 
for  well  over  a  decade. 
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The  detection  of  G1  876d  by  Rivera  et  al.  (2005)  underscored 
the  rapid  development  of  the  Doppler  technique.  The  detection 
of  planets  having  only  a  few  Earth  masses  is  the  culmination  of 
advances  described  by  Butler  et  al.  (1996),  Marcy  &  Butler  (1998), 
Marcy  et  al.  (2000,  2005),  and  Lovis  et  al.  (2006).  Neverthe¬ 
less,  while  radial  velocity  precision  has  improved  to  better  than 
1  m  s^',  Keplerian  orbital  fits  to  Doppler  velocities  give  estimates 
of  A/pisin  (i)  rather  than  Mpi.  In  order  to  ascertain  the  true  mass 
of  a  planet,  one  must  obtain  an  independent  measure  of  the  orbital 
inclination  i,  which  is  most  easily  measured  when  the  planet  is 
observed  to  transit  across  the  face  of  the  parent  star. 

The  now  celebrated  companion  to  HD  209458  provided  the 
first  example  of  a  transiting  extrasolar  planet,  and  this  body  has 
generated  an  incredible  variety  of  observational  and  theoretical  in¬ 
vestigations.  HD  209458b  was  initially  detected  with  the  radial 
velocity  method,  and  was  then  discovered  to  transit  on  the  basis 
of  follow-up  photometric  observations  at  the  predicted  transit  times 
(Henry  et  al.  1999, 2000;  Charbonneau  et  al.  2000).  HD  209458b ’s 
well-sampled  light  curve,  when  combined  with  high-precision 
radial  velocity  data,  has  permitted  precise  measurements  of  the 
planet’s  basic  attributes  (e.g.,  Mazeh  et  al.  2000;  Brown  et  al. 
2001;  Laughlin  et  al.  2005;  Winn  &  Holman  2005). 

To  date,  nine  additional  planets  have  been  found  to  transit  their 
parent  stars,  and  the  varied  uses  to  which  these  transits  have  been 
put  are  described  in  the  review  article  by  Charbonneau  et  al.  (2006). 
Our  goal  in  this  paper  is  to  ascertain  whether  the  two  outer  planets 
orbiting  G1  876  (IL  Aqr,  GJ  876,  HIP  113020)  can  similarly  be 
observed  in  transit.  A  positive  detection  would  represent  a  further 
major  improvement  in  our  understanding  both  of  this  specific 
system  and  of  planetary  properties  in  general. 

Rivera  et  al.  (2005)  carried  out  a  photometric  search  for  tran¬ 
sits  by  planet  d  (P  ^  1.974  day),  and  showed  that  such  transits 
do  not  occur.  To  date,  photometric  searches  for  transits  by  planet 
c  (P  ^  30  day)  and  b  (P  ~  60  day)  have  not  been  reported.  Al¬ 
though  the  a  priori  geometric  transit  probabilities  for  G1 876c  and 
b  are  only  Ptr  =  1 .4%  and  Ptr  =  0.9%,  respectively,  the  expected 
photometric  transit  depths  are  greater  than  10%.  Librational  mo¬ 
tions  arising  from  the  resonance  between  the  two  outer  planets 
lead  to  predicted  transit  epochs  for  the  outer  planets  that  are  not 
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TABLE  1 

CoPLANAR  Three-Planet  i  —  90°  Jacobi  Parameters  of  the  G1  876  System 


Parameter 

Planet  d 

Planet  c 

Planet  b 

. 

5.89  ±  0.54  Mg, 

0.619  ±  0.005Mjup 

1.935  ±  0.007Mjup 

P  (days) . 

1.93776  ±  0.00007 

30.340  ±  0.013 

60.940  ±  0.013 

K  (m  s-') . 

6.46  ±  0.59 

88.36  ±  0.72 

212.60  ±  0.76 

riGAU) . 

0.0208067  ±  0.0000005 

0.13030  ±  0.00004 

0.20783  ±  0.00003 

e . 

0  (fixed) 

0.2243  ±  0.0013 

0.0249  ±  0.0026 

at  (deg) . 

0  (fixed) 

198.3  ±  0.9 

175.7  ±  6.0 

M  (deg) . 

309.5  ±  5.1 

308.5  ±  1.4 

175.5  ±  6.0 

transit  epoch . 

JD  2,452,490.756  ±  0.027 

JD  2,452,517.633  ±  0.051 

^  Quoted  uncertainties  in  planetary  masses  and  semimajor  axes  do  not  incorporate  the  uncertainty  in  the  mass  of  the 
star. 


spaced  evenly  in  time.  (The  resonant  structure  of  the  system  is  dis¬ 
cussed  further  in  Laughlin  &  Chambers  200 1 ;  Rivera  &  Lissauer 
2001;  Kley  2003;  and  Beauge  et  al.  2005). 

Benedict  et  al.  (2002)  used  Hubble  Space  Telescope  (HST) 
astrometry  to  derive  an  orbital  inclination  for  G1  876b  of  i  = 
84°  ±  6° .  Their  result  thus  suggests  that  planet  b  has  a  mass  equal 
or  close  to  the  radial  velocity  ( RV )  -  determined  value  of  M  sin  (i)  = 
1.89Mjup,  and  that  the  system  is  within  1  ct  of  being  viewed 
edge-on.  The  implication  of  the  Benedict  et  al.  (2002)  mea¬ 
surement  is  that  occurrence  of  transits  is  much  more  likely  than 
the  a  priori  geometric  probability  would  suggest.  The  Benedict 
et  al.  (2002)  result  is  in  conflict,  however,  with  the  dynamical 
models  of  the  system  presented  by  Rivera  et  al.  (2005),  who  And 
a  best  fit  to  the  radial  velocity  data  when  the  system  (assumed  co- 
planar)  is  inclined  by  i  ~  50°  ±  3°  with  respect  to  the  plane  of  the 
sky.  Therefore,  if  transits  are  detected,  then  the  Benedict  et  al. 
(2002)  results  will  be  confirmed  at  the  expense  of  the  conclu¬ 
sions  of  Rivera  et  al.  (2005). 

G1 876b  and  G1 876c  are  members  of  an  as  yet  unobserved  class 
of  giant  planets  with  intermediate  surface  temperatures.  The  ex¬ 
pected  effective  temperatures  of  7);^  b  =  160  K  and  7).^  c  =  200  K 
of  the  planets  are  much  lower  than  the  value  ~  1350  K  mea¬ 
sured  with  Spitzer  for  HD  209458b  (Charbonneau  et  al.  2005), 
but  still  considerably  higher  than  the  T’^jf  ^  135  K  measured  for 
Jupiter.  Determination  of  their  physical  properties  would  thus  pro¬ 
vide  a  useful  link  between  Jupiter  and  Saturn  on  the  one  hand,  and 
objects  such  as  HD  209458b  on  the  other.  Accurate  measurement 
of  their  radii  would  also  give  clues  to  their  interior  structure  (see 
Charbonneau  et  al.  2006). 

Finally,  the  G1  876  system  with  its  strongly  gravitationally  in¬ 
teracting  planets  and  large  planet-to-star  radius  ratios,  has  the 
potential  to  display  an  extraordinarily  informative  set  of  transit 
light  curves.  Our  model  light  curves  here  show  what  one  can  ex¬ 
pect  when  interacting  transiting  systems  are  discovered,  for  ex¬ 
ample,  with  ongoing  radial  velocity  surveys  or  with  space-based 
missions  such  as  Kepler  (Borucki  et  al.  2003)  or  COROT  (Baglin 
2003),  and  allow  for  a  more  complete  understanding  of  the  multi¬ 
planet  dynamics. 

The  plan  of  this  paper  is  as  follows.  In  §  2  we  describe  our  dy¬ 
namical  model  of  fhe  system,  and  illustrate  some  of  the  photo¬ 
metric  properties  that  the  in-transit  light  curves  would  be  expected 
to  display.  In  §  3  we  describe  an  evaluation  of  the  Doppler  radial 
velocity  measurements  that  have  been  taken  during  predicted  tran¬ 
sit  intervals.  In  §  4  we  evaluate  a  range  of  archival  and  newly 
obtained  photometric  data  sets  for  evidence  of  transits.  In  §  5  we 
discuss  our  results,  as  well  as  the  current  status  and  outlook  for 
the  ongoing  Transitsearch  collaboration  and  conclude  in  §  6. 


2.  THE  DYNAMICAL  MODEL 

G1 8 76  ( M4 V )  is  the  40th  nearest  star  to  Earth,  with  a  Hipparcos- 
determined  distance  of  4.69  pc  (Perryman  et  al.  1997);  Tycho-II 
and  UCAC  positions  validate  this  determination.  As  discussed  in 
Laughlin  et  al.  (2005),  current  estimates  of  its  mass  and  radius 
stand  at  M*  =  0.32  Mq  and  R*  =  0.3  TJq.  We  adopt  these  val¬ 
ues  as  fixed  in  this  paper.  To  date,  1 55  Doppler  velocities  mea¬ 
sured  at  the  Keck  telescope  and  16  velocities  measured  at  the 
Lick  telescope  have  been  published  for  the  G1 876  system  (Marcy 
et  al.  200 1 ;  Rivera  et  al.  2005).  The  two  planets  induce  a  total  ve¬ 
locity  half-amplitude  of  the  star  of  nearly  0.25  km  s^*.  Further¬ 
more,  as  discussed  in  Rivera  et  al.  (2005),  the  system  displays  a 
very  low  level  of  stellar  radial  velocity  “jitter,  ”  which  further  aids 
in  obtaining  a  detailed  characterization  of  the  orbits.  The  system 
has  now  been  observed  for  more  than  80  orbital  periods  of  the 
middle  planet,  c.  This  extensive  time  baseline  allows  for  a  much 
more  detailed  study  of  planet-planet  interactions  than  can  be  ob¬ 
tained  with  the  other  exoplanetary  systems  known  to  be  in  2 : 1 
resonance,  a  list  that  now  includes  HD  82943  (Mayor  et  al.  2004; 
Lee  et  al.  2006),  HD  128311  (Vogt  et  al.  2005),  and  HD  73526 
(Tinney  etal.  2006).  These  three  systems  all  have  inner  planet  pe¬ 
riods  on  the  order  of  10  times  longer  than  G1  876c.  Our  baseline 
coplanar  i  —  90°  orbital  model  was  obtained  by  Rivera  et  al.  (2005) 
from  a  dynamical  fit  to  the  Keck  radial  velocity  data.  For  reference, 
the  parameters  of  this  model  are  reproduced  in  Table  1 .  Additional 
fits  with  i  <  90°  were  also  obtained  to  the  Keck  data  using  the  al¬ 
gorithm  described  by  Rivera  et  al.  (2005).  As  with  the  i  =  90° 
baseline  model,  the  system  is  assumed  to  contain  three  planets, 
all  orbiting  the  star  in  the  same  plane.  We  varied  the  inclination 
of  the  orbital  plane  to  the  plane  of  the  sky  from  i  =  90°  to  89?0  in 
decrements  of  0?1,  and  fitted  for  each  of  the  13-1-1  remaining 
parameters.  Then  V-body  integrations  of  the  fit  generate  a  sequence 
of  transit  predictions  for  the  outer  two  planets.  Uncertainties  in 
the  predicted  transit  midpoints  were  generated  with  the  bootstrap 
Monte  Carlo  technique  described  by  Rivera  et  al.  (2005).  For  both 
outer  planets,  the  uncertainty  in  the  transit  midpoint  times  was 
found  to  vary  from  ~1  to  a  few  hours.  The  variation  in  uncer¬ 
tainty  from  transit  to  transit  is  plotted  in  Figure  1.  The  predicted 
transit  midpoints  are  best  constrained  during  the  epoch  spanned 
by  the  radial  velocity  observations.  Errors  increase  as  the  transit 
predictions  are  generated  for  earlier  or  later  times.  Light  curves 
were  produced  with  a  photometric  model  of  the  transit,  which  as¬ 
sumes  that  the  planets  are  opaque  disks  with  radii  R  =  0.937?jup 

Note  that  the  elements  reported  in  Table  1  are  osculating  orbital  elements 
(see,  e.g.,  Murray  &  Dermott  1999)  referenced  to  the  epoch  JD  2,452,490. 
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Fig.  1. — Uncertainty  in  the  central  transit  time  for  planets  G1  876b  {solid  line) 
and  G1 876c  {dashed  line).  The  vertical  dashed  lines  delineate  the  time  span  of  the 
Keck  RV  observations  used  to  generate  the  coplanar  fits  that  form  the  basis  of  our 
model  photometric  and  RV  curves.  The  vertical  dotted  line  indicates  the  epoch  for 
the  fits  (JD  2,452,490.0). 

for  c  and  7?  =  1 .047?jup  for  b.  These  equatorial  radii  are  computed 
using  the  models  of  Bodenheimer  et  al.  (2003).  The  specific  in¬ 
tensity  of  the  stellar  disk  is  modeled  with  a  linear  limb  darkening 
law,  given  as 

=  (1) 

where  I  is  the  specific  intensity  of  the  stellar  disk,  v=  0.724 
is  the  corresponding  limb  darkening  coefficient  appropriate  to 
K-band  observations  of  an  M4V  primary  (per  Claret  2000),  and 
^  =  cos  (9). 

Figures  2  and  3  show  a  selection  of  model  light  curves  for 
planets  b  and  c.  For  planet  b,  which  has  an  eccentric  (e  =  0.22), 
rapidly  precessing  orbit  (dzuldt  ^  44°  yr^^),  the  shape  of  the 


t  (JD  -  2450000) 

Fig.  2. — Model  limb-darkened  light  curves  during  a  predicted  transit  by 
planet  G1 876b  when  it  is  near  periastron  {solid  curves)  and  near  apastron  {dotted 
curves).  The  curves  correspond  to  different  assumed  inclinations  {i  —  90°,  89?9, 
89?8,. . .,  89?5).  The  shallowest  occur  for  i  —  89?5. 


t  (JD  -  2450000) 

Fig.  3 . — Model  transit  light  curves  for  planet  G1 876c,  when  it  is  near  periastron 
{solid  curves),  and  near  apastron  {dotted  curves).  Inclinations,  i  —  90°  through 
i  —  89?0  are  shown  in  decrements  of  0?1.  Planet  G1  876c ’s  eccentricity  e  ~  0.22 
leads  to  a  greater  difference  in  transit  depth  and  duration  between  periastron  and 
apastron. 

predicted  curves  depends  strongly  on  the  osculating  periastron 
angle.  For  i  <  90°,  the  transits  are  both  shorter  and  deeper  when 
the  transits  occur  near  periastron.  The  orbital  eccentricity  of  planet 
b  is  more  nearly  circular,  with  e  =  0.025.  In  this  case,  there  is  little 
variation  in  the  predicted  light  curves  as  the  orbits  process.  The 
known  transiting  planets  have  transit-to-transit  intervals  AT)  that 
are  evenly  spaced  to  within  the  current  resolution  of  the  obser¬ 
vations.  HD  209458b,  for  example,  has  had  individual  transit  mid¬ 
points  measured  to  an  accuracy  of  several  seconds  (Brown  et  al. 
2001;  Wittenmyer  et  al.  2005),  and  an  average  period  measured 
to  an  accuracy  of  83  ms  (Wittenmyer  et  al.  2005;  Winn  &  Holman 
2005).  If,  however,  one  is  able  to  measure  differences  in  the  time 
intervals  between  successive  transits,  one  can  derive  additional 
dynamical  information.  As  discussed  by  Miralda-Escude  (2002), 
Holman  &  Murray  (2005),  and  Agol  et  al.  (2005),  measurements 
of  small  variations  from  exact  periodicity  can  be  used  to  infer  the 
presence  of  additional  bodies  in  the  system,  potentially  with  the 
masses  in  the  range  appropriate  to  terrestrial  planets.  In  the  event 
that  either  G1 876b  or  G1 876c  were  to  transit,  one  would  expect 
the  departures  from  strict  periodicity  to  be  large.  In  fact,  succes¬ 
sive  transits  of  planet  c  would  occur  at  intervals  differing  by  as 
much  as  4  hr.  By  fitting  a  dynamical  model  to  such  transit  inter¬ 
vals,  one  could  obtain  information  regarding  all  of  the  orbital 
elements  of  the  planets  (including  inclinations  and  nodes). 

3.  SEEKING  EVIDENCE  OF  TRANSITS  IN  THE  RADIAL 
VELOCITY  DATA  SET 

During  a  planetary  transit,  the  occulting  disk  of  the  planet 
passes  over  portions  of  the  stellar  disk  containing  a  varying 
component  of  line-of-sight  velocity  arising  from  stellar  rotation; 
this  phenomenon  produces  an  asymmetry  in  the  star’s  spectral 
lines  known  as  the  Rossiter-McLaughlin  effect  (Rossiter  1924; 
McLaughlin  1924).  The  effect  has  been  described  in  detail  in  the 
context  of  planetary  transits  by  Ohta  et  al.  (2005),  and  has  been 
analyzed  for  HD  209458  by  Queloz  et  al.  (2000),  Bundy  &  Marcy 
(2000),  and  Winn  et  al.  (2005),  and  for  HD  149026  by  Wolf  et  al. 
(2006).  In  fact,  the  transiting  planet  HD  1 89733  was  initially  dis¬ 
covered  to  transit  not  via  photometry,  but  rather  via  such  in-transit 
radial  velocity  variations  found  by  Bouchy  et  al.  (2005).  For  our 
study  of  G1  876,  we  compared  the  observation  epochs  of  all  of 
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Fig.  4. — Three-planet  RV  model  (dashed  line)  with  a  superimposed  model 
Rossiter  curve  (solid  line)  during  a  predicted  transit  by  planet  G1 876c  near  epoch 
JD  2,453,301.79.  The  observed  Keck  RV  data  are  shown  as  the  black  filled-in 
circles  with  error  bars. 

the  Lick  (see  Marcy  et  al.  2001)  and  Keck  (see  Rivera  et  al.  2005) 
radial  velocity  observations  to  see  if  any  of  these  observations  oc¬ 
curred  during  a  predicted  transit  (assuming  i  =  90°).  For  G1 876, 
we  used  a  rotational  period  that  was  determined  to  be  96.7  days 
using  rotational  modulation  (as  reported  in  Rivera  et  al.  2005); 
we  used  this  value  to  calculate  the  amplitude  of  the  Rossitter- 
McLaughlin  effect.  Four  such  observations  were  found,  three  taken 
during  a  predicted  transit  by  planet  c,  and  one  taken  during  a  pre¬ 
dicted  transit  by  planet  b.  F igures  4  and  5  display  the  relevant  parts 
of  the  radial  velocity  model  (both  with  and  without  the  Rossiter- 
McLaughlin  effect),  with  the  observations  shown.  The  radial  ve¬ 
locity  model  curves  in  both  figures  were  generated  by  adding 
template  curves  that  use  the  stellar  and  planetary  radii  and  limb 
darkening  law  given  above  (and  that  apply  the  analytic  integrals 
given  in  Ohta  et  al.  2005),  to  our  fiducial  i  =  90°  three-planet 
radial  velocity  model.  For  the  transit  epoch  in  the  case  of  planet  c 
(Fig.  4),  the  data  are  inconsistent  with  a  transit.  In  Figure  5,  for  the 
outer  planet  b,  the  situation  is  unresolved. 


Fig.  6. — Model  i  —  89?6  light  curves  for  G1  876b  and  c  superimposed  on 
Hipparcos  epoch  photometry.  None  of  the  Hipparcos  measurements  were  made 
near  or  during  a  predicted  transit  window. 

4.  A  SEARCH  FOR  TRANSITS  IN  PHOTOMETRIC 
DATA  SETS 

We  have  also  compared  our  model  light  curves  to  a  variety  of 
photometric  data  sets  for  G1 876  taken  since  1989.  The  results  of 
these  comparisons  are  as  follows. 

Hipparcos. —  The  Hipparcos  mission  (see  Perryman  et  al.  1 997) 
produced  67  accepted  photometric  measurements  for  G1  876 
(HIP  113020).  These  measurements  have  a  median  magnitude 
Hp  =  10.148,  with  (7  =  0. 0045 .  There  is  evidence  for  a  weak  peri¬ 
odicity  when  the  data  is  folded  at  the  96.7  days  rotational  period  of 
the  star  (also  noted  in  Rivera  et  al.  2005).  In  Figure  6,  we  plot  our 
model  i  =  89f6  light  curves  for  G1  876b  and  c  superimposed  on 
the  Hipparcos  epoch  photometry.  None  of  the  Hipparcos  mea¬ 
surements  were  made  near  or  during  a  predicted  transit  window. 

Fairborn  Observatory . —  A  number  of  long-term  photometric 
monitoring  programs  are  being  carried  out  with  the  automated 
telescopes  at  Fairborn  Observatory  (e.g.,  Henry  1999;  Eaton  et  al. 
2003)  and,  as  described  in  detail  in  Rivera  et  al.  (2005),  G1 876  has 
been  on  one  of  the  observing  programs.  Figures  7  and  8  show  the 
Fairborn  Observatory  data  and  our  corresponding  model  light 


t  (JD  -  2452446.0) 


Fig.  5. — Three-planet  RV  model  {dashed  line)  with  a  superimposed  model 
Rossiter  curve  {solid  line)  during  a  predicted  transit  by  planet  G1 876b  near  epoch 
JD  2,452,446.14.  Again,  the  observed  Keck  RV  data  are  shown  as  the  black  filled- 
in  circle  with  error  bars. 


Fig.  7. — Model  i  —  89?3  light  curves  for  G1  876b  and  c  superimposed  on 
Fairborn  Observatory  photometry.  Note  that  the  intrinsic  long-term  variability  for 
G1  876  over  a  ~  1000  day  period  is  of  order  0.05  mag. 
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Fig.  8. — Model  light  curves  for  G1  876c  superimposed  on  Fairborn  Obser¬ 
vatory  photometry  taken  near  three  individual  predicted  transit  windows.  The 
solid  lines  correspond  to  model  light  curves  arising  from  assumed  inclinations 
running  from  i  —  90?0  to  i  —  89?0.  The  light  lines  correspond  to  i  —  90°  model 
light  curves  arising  from  bootstrap  fits  in  which  the  transit  is  3  cr  early,  1  u  early, 
and  1  (7  and  3  u  late. 


curve  for  planets  b  and  c.  Data  taken  during  predicted  transit  epochs 
for  planet  c  show  no  indication  of  a  transit.  The  Fairborn  Obser¬ 
vatory  data  indicate  that  long-term  stellar  variability  from  G1 876 
is  of  the  order  0.05  mag,  which  is  somewhat  smaller  than  the 
variation  observed  with  Hipparcos. 

Transitsearch  and  AAVSO^^. —  Seagroves  et  al.  (2003)  de¬ 
scribe  Transitsearch  as  a  cooperative  distributed  observing 
project  involving  submeter  class  telescopes  worldwide.  The  Tran¬ 
sitsearch  strategy  is  to  observe  known  planet-bearing  stars  at  the 
dates  and  times  when  transits  are  expected  to  occur.  It  is  therefore 
a  targeted  search,  which  differentiates  it  from  ongoing  wide-held 
surveys.  We  identify  observing  windows  for  candidate  stars 
through  the  use  of  the  bootstrap  Monte  Carlo  technique  described 
by  Laughlin  et  al.  (2005).  The  observational  campaigns  are  pri- 


See  http://www.aavso.org. 


oritized  by  the  a  priori  likelihood  that  a  particular  candidate  planet 
will  display  transits.  This  likelihood  is  given  by 

P,„,  =  0.0045  (ift)  (^4^)  , 

(2) 

where  a  is  the  semimajor  axis  of  the  orbit,  7?*  is  the  radius  of 
the  star,  R^i  is  the  radius  of  the  planet,  va  is  the  longitude  of 
periastron,  and  e  is  the  orbital  eccentricity. 

The  hrst  Transitsearch  data  set  for  G1 876  was  obtained  by  one 
of  us  (P.  D.  S.  Shankland)  from  2003  October  27.9583  to  October 
28.0368,  with  a  Dhltered,  16-bit  speckle  CCD  at  the  primary  focus 
of  a  40.4  cm  f/5.1  rehector  (at  air  mass  from  secz  =  2. 7-1. 7). 
The  photometric  time  series  shows  a  monotonic  hux  increase  of 
0.16  ±  0.04  mag  from  2003  October  27.9583  to  October  27.9972. 
As  shown  on  the  left  panel  of  Figure  9,  the  depth,  duration,  and 
timing  of  this  event  were  consistent  with  an  egress  from  transit  of 
G1 876c,  in  agreement  with  our  predictions.  The  possibility  that  a 
transit  had  been  observed  spurred  us  to  orchestrate  a  distributed 
observational  campaign  during  the  2004  season. 

During  the  2004  June  24  opportunity,  observations  were  at¬ 
tempted  by  10  observers  in  Australia.  Poor  weather  thwarted  all 
but  one  observer,  and  no  transit  signal  was  seen  in  the  data,  which 
covered  about  half  of  the  1  a  transit  window.  Eight  Australian  and 
Japanese  observers  participated  in  the  following  2004  July  24 
planet  c  campaign.  They  were  all  clouded  out.  Weather  improved 
for  the  2004  August  23  planet  c  opportunity,  permitting  obser¬ 
vations  by  five  of  the  nine  Australian,  South  African,  and  German 
observers.  No  transit  signal  was  seen  with  photometric  depth  greater 
than  1%. 

In  2004  October,  a  number  of  AAVSO  and  Transitsearch  ob¬ 
servers  obtained  photometry  during  a  closely  spaced  pair  of  tran¬ 
sit  windows  for  planets  b  and  c.  Begiiming  October  20  1 1 : 22  (UT) 
and  ending  October  23  15:27  (UT),  a  total  of  2795  CCD  ob¬ 
servations  were  obtained.  In  addition,  2981  photometric  observa¬ 
tions  were  obtained  in  the  days  immediately  before  and  after  the 
window  to  set  baseline  activity  and  to  look  for  red  dwarf  flaring. 
Such  a  large,  distributed  network  of  observers  allowed  coverage 
of  the  window  with  very  few  observing  gaps  over  the  3  day  period. 
Gaps  in  the  coverage  were  smaller  than  the  estimated  transit  period 
or  forecast  window,  allowing  us  to  conclude  that  transits  for  planets 
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Fig.  9. — Left:  data  for  G1876c,  this  plate  shows  dynamical  model  light  curves  again  for  the  same  decrements  of  i,  taken  near  JD  2,452,940.4381;  relative  photometry, 
which  began  in  asfronomical  twilight  {filled  circles).  Right:  Dynamical  model  light  curves  for  G1 876  due  to  planet  c  for  i  —  90?0  through  89?0  in  0?1  decrements,  taken  near  a 
predicted  transit  centered  at  JD  2,453,301.7879.  The  data  were  compiled,  baselined,  and  shifted  to  employ  a  mean  magnitude  of  F  =  10.19,  based  on  relative  photometry 
obtained  by  AAVSO  and  Transitsearch  {filled  circles).  The  light  lines  correspond  to  i  —  90°  model  light  curves  arising  from  bootstrap  fits  in  which  the  transit  is  3  cr  early,  1  cr 
early,  and  1  cr  late  and  3  cr  late. 
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Fig.  10. — Dynamical  model  light  curves  for  G1  876  due  to  planet  b  for  i  — 
90?0  through  89?5  in  0?1  decrements,  near  a  predicted  transit  centered  at  JD 
2,453,300.3806.  As  with  the  previous  hgure,  the  data  were  baselined  and  shifted 
to  employ  a  mean  magnitude  of  V  —  10.19,  based  on  relative  photometry  ob¬ 
tained  by  AAVSO  and  Transitsearch  {filled  circles).  The  light  lines  correspond 
to  i  —  90°  model  light  curves  arising  from  bootstrap  fits  in  which  the  transit  is 
3  cr  early,  1  u  early,  and  1  cr  late  and  3  cr  late. 

b  and  c  did  not  occur  during  the  2004  October  campaign.  Photo¬ 
metric  data  from  the  campaign  are  plotted  in  Figures  9  and  10. 

Individual  observations  and  uncertainty  estimates  from  the 
2004  October  campaign  are  available  from  the  AAVSO  Web  site, 
while  the  other  data  sets  can  be  found  at  the  Transitsearch 
website.'^  We  are  continuing  to  obtain  photometry  of  G1  876, 
with  the  goal  of  gaining  a  better  understanding  of  the  photometric 
variability  of  mature  (age  exceeding  1  Gyr)  red  dwarfs. 

5.  DISCUSSION 

The  Transitsearch  photometric  network  has  been  in  operation 
for  3  yr,  and  has  participating  observers  capable  of  providing 
fully  global  longitude  and  latitude  coverage.  The  network  main¬ 
tains  a  continuously  updated  catalog'®  of  the  known  census  of  ex¬ 
trasolar  planets.  The  catalog  is  currently  the  only  available  source 
for  (1)  a  priori  transit  probabilities  and  predicted  transit  depths  for 
the  known  extrasolar  planets,  (2)  predicted  transit  ephemeris  and 
transit  window  uncertainties  based  on  orbital  fits  to  published  radial 
velocities,  and  (3)  the  results  of  known  photometric  searches  for 
transits  of  planet-bearing  stars.  To  date,  the  majority  of  negative 
transit  searches  have  not  been  reported  in  the  literature.  With  this 
paper,  we  hope  to  spur  a  reversal  of  that  trend. 

Due  to  significant  uncertainties  in  the  orbital  fits,  it  is  generally 
very  dilficult  to  definitively  rule  out  transits  by  planets  with  P  > 
10  day.  The  G1  876  planets  form  a  notable  exception  to  this  rule 
because  of  their  high  signal-to-noise  ratio  T  ~  50,  and  the  exten¬ 
sive  radial  velocity  data  set  allows  for  a  very  accurate  orbital  fit. 
Campaigns  by  the  Transitsearch  collaboration  have  thus  far  led 
to  transits  of  a  number  of  other  shorter-period  planet-bearing  stars 
being  ruled  out  to  a  high  degree  of  confidence.  The  planets  for 
which  the  possibility  of  transits  have  been  significantly  discounted 
by  our  campaigns  include  FID  217107b  (P  =  7.127  day,  Ptr  = 
8.0%;  see  also  Vogt  et  al.  2005),  HD  168746b  {P  =  6.403  day, 


Our  AAVSO  data  are  available  at  http://www.aavso.org/data/download/, 
while  other  collected  data  can  be  found  at  http://www.ucolick.org/%7elaugh/ 
G1876/ _ c.results.html. 

Available  at  http://www.ucolick.org/%7elaugh/. 


=  8.1%;  see  also  Pepe  et  al.  2002),  and  HD  68988b  {P  = 
6.276  day,  P^^  —  9.0%;  see  also  Vogt  et  al.  2002). 

Furthermore,  planets  for  which  usable  photometry  of  the  parent 
stars  has  been  obtained  by  our  network  during  the  3  a  transit  win¬ 
dows,  but  for  which  transits  cannot  yet  be  fully  ruled  out,  include 
HD  188753 Ah  {P  =  3.348  day,  Htr  =  1 1  -8%),  HD  76700b  {P  = 
3.971  day,  P^  =  10.0%),  HD  13445b  {P  =  15.76  day,  = 
3.4%),  HD  74156b  (P  =  51.64  day,  P^  =  3.8%),  HD  37605b 
{P  =  54.23  day,  Ptr  =  2.2%h  and  HD  80606b  (P  =  1 1 1.4  day, 
Ptr  =  1.7%).  Planets  for  which  campaigns  are  scheduled  to  be¬ 
gin  soon  include  GL  581b  (P  =  5.366  day, Ptr  =  3.6%),andHD 
99492b  (P  =  17.04  day,  =  3.4%).  GL  581  is  of  particular 
interest  because  the  primary  star  is  an  M  dwarf,  and  the  planet  (if 
it  transits)  will  be  a  Neptune-mass  object  M  sin  i  —  0.052Mjup 
(Bonfils  et  al.  2005).  The  1.6%  transit  depth  would  allow  an  un¬ 
precedented  physical  characterization  of  a  low-mass  extrasolar 
planet.  To  date,  our  G1 876  campaign  has  employed  the  most  glob¬ 
ally  far-reaching  distribution,  and  along  with  recent  leaps  in  fidel¬ 
ity  in  similar  campaigns,  clearly  underscores  the  utility  of  the 
network. 

What  is  more,  G1 876b  and  c  constitute  a  dramatic  exception  to 
the  emerging  aphorism  that  Jupiter-mass  planets  are  rarely  asso¬ 
ciated  with  M  dwarfs.  Indeed,  a  paucity  of  giant  planets  orbiting 
red  dwarfs  seems  to  be  a  natural  consequence  of  the  core  accre¬ 
tion  theory  of  planet  formation  (Laughlin  et  al.  2004).  The  core¬ 
accretion  theory  predicts,  however,  that  Neptune-mass  and  smaller 
mass  planets  will  be  very  common  around  red  dwarfs.  This  stands 
in  contrast  to  the  predictions  of  Boss  (2006),  who  suggests  that 
the  gravitational  instability  mechanism  is  the  dominant  mode  of 
giant  planet  formation.  If  this  is  true,  then  Jupiter-mass  planets 
should  be  just  as  common  around  red  dwarfs  as  they  are  around 
solar-type  stars.  Efforts  to  detect  transits  of  small  planets  orbiting 
small  stars  are  therefore  likely  to  gain  increasing  importance  and 
focus  over  the  next  several  years. 

There  are,  in  fact,  several  observational  indications  that  low- 
mass  planets  may  be  common  around  red  dwarfs.  Radial  velocity 
surveys  have  recently  reported  the  detection  of  Neptune-mass  com¬ 
panions  in  short-period  orbits  around  the  red  dwarfs  GL  581  and 
GL  436.  The  Optical  Gravitational  Lensing  Experiment  (OGLE) 
team  of  observers  using  the  microlensing  method  have  detected 
the  signature  of  what  appears  to  be  a  planet  with  5.5  Mq  orbiting 
a  distant  red  dwarf  (Beaulieu  et  al.  2006).  With  these  points  in 
mind,  we  can  put  the  results  of  our  G1  876  campaign  in  a  broader 
context.  A  ground-based  (and  “fiscally  viable”)  photometric  net¬ 
work  such  as  Transitsearch,  especially  in  collaboration  with  a  de¬ 
dicated  bank  of  telescopes,  can  capably  monitor  individual 
M  dwarfs  to  search  for  transits  of  terrestrial-sized  planets.  Owing 
to  the  intrinsic  long-term  variability  of  M  dwarfs,  candidate  plan¬ 
ets  would  have  to  be  identified  on  the  basis  of  a  full  transit  sig¬ 
nature  in  a  single  time  series. 

A  0.1  Mq  M  dwarf  has  RJRq^QA,  T^s^llSQ)  K,  and 
LJLq^  5  X  10^^;  a  habitable  planet  that  receives  an  Earth- 
equivalent  fiux  from  the  star  thus  needs  to  orbit  at  a  distance  of 
just  0.022  AU,  which  corresponds  to  an  orbital  period  of  3.85  days. 
An  Earth-like  planet  with  such  a  period  would  be  rotationally 
synchronized  to  the  red  dwarf,  and  in  the  absence  of  any  signif¬ 
icant  perturbing  bodies,  its  orbit  would  be  almost  perfectly  cir¬ 
cular.  Interestingly,  however,  simulations  show  that  habitability 
is  unlikely  to  be  adversely  affected  by  a  spin-orbit  period  syn¬ 
chronization.  Joshi  et  al.  (1997)  used  a  global  climate  model  to 
investigate  how  the  Earth’s  climate  would  respond  if  the  Earth 
were  tidally  locked  to  the  Sim,  and  found  that  Earth  remains  habit¬ 
able  in  this  configuration,  at  1  AU,  and  perhaps  throughout  the 
habitable  zone. 
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TheM  =  7.5  M^,P  =  1.9379  day  companion  to  G1 876  dem¬ 
onstrates  that  it  is  not  unreasonable  to  expect  terrestrial-mass  bodies 
on  short-period  orbits  of  0.1  Mq  stars.  R.  Montgomery  &  G. 
Laughlin  (private  communication,  2006)  have  carried  out  accre¬ 
tion  simulations  with  the  Wetherill-Chambers  method  (  Wetherill 
1996;  Chambers  2001)thatmodel  the  accretion  of  terrestrial-mass 
planets  in  short-period  orbits  about  M  dwarfs.  These  simulations 
provide  further  evidence  that  terrestrial-mass  planets  will  commonly 
form  in  the  desired  short-period  orbits. 

If  habitable  planets  do  commonly  form  in  orbit  around  low 
mass  M  dwarfs,  the  chances  of  detecting  and  characterizing  them 
are  surprisingly  good.  The  transit  of  an  Earth-sized  planet  will 
block  about  1%  of  the  stellar  flux  of  a  0.1  A/®  star.  For  a  planet 
on  a  habitable  3.85  day  orbit,  the  transit  will  be  relatively  brief, 
^40  minutes.  The  a  priori  geometric  probability  of  observing 
such  a  transit  is  2%. 

A  1%  photometric  dip  is  readily  detectable.  Amateur  astrono¬ 
mers  who  participate  in  the  Transitsearch  collaboration  routinely 
achieve  detection  thresholds  of  considerably  better  than  1  %,  as 
evidenced  by  confirming  detections  of  HD  149026b,  which  has 
a  transit  depth  of  just  0.3%.  Indeed,  capable  amateur  observers 
have  demonstrated  the  photometric  capability  to  detect  the  pas¬ 
sage  of  a  Mars-sized  body  in  front  of  an  11th  mag  0.1  Mq  red 
dwarf 

Inevitably  one  must  ask  how  many  suitable  red  dwarfs  are 
available  on  the  sky  as  a  whole.  Even  though  the  lowest  mass 
M  stars  are  the  most  common  type  of  star,  they  are  also  exceedingly 
dim.  A  dedicated  3  m  class  telescope  would  be  required  to  prop¬ 
erly  search  the  magnitude  range  V  =  16-18,  where  observations 
of  the  lowest  mass  red  dwarfs  begin  to  be  plentiful  (and  indeed, 
many  0.1  Mq  stars  within  10  pc  still  likely  remain  to  be  discov¬ 
ered).  In  the  meantime,  however,  we  recommend  that  observers 
obtain  high-cadence  photometry  of  nearby  single  stars  listed  in 
the  RECONS  catalog'^  of  the  100  nearest  stellar  systems;  these 
include  Proxima  Centauri,  Barnard’s  Star,  Wolf  359,  Ross  154, 
Ross  128,  DX  Cancri,  GJ  1061,  GJ  54.1,  and  GJ  83.1.  We  esti¬ 
mate  that  each  one  has  a  ~  1  %  chance  of  harboring  a  detectable 
transiting,  potentially  habitable  planet,  readily  detectable  by  dis¬ 
tributed  photometric  observation. 


”  Maintained  by  T.  Henry  and  collaborators  at  Georgia  State  University. 


Agol,  E.,  Steflfen,  J.,  Sari,  R.,  &  Clarkson,  W.  2005,  MNRAS,  359,  567 
Baglin,  A.  2003,  Adv.  Space  Res.,  31,  345 

Beauge,  C.,  Michtchenko,  T.,  &  Ferraz-Mello,  S.  2006,  MNRAS,  365,  1160 

Beaulieu,  J.-R,  et  al.  2006,  Nature,  439,  437 

Benedict,  G.  R,  et  al.  2002,  Apt,  581,  LI  15 

Bodenheimer,  R,  Laughlin,  G.,  &  Lin,  D.  N.  C.  2003,  Apt,  592,  555 

Bonfils,  X.,  et  al.  2005,  A&A,  443,  L15 

Borucki,  W.  J.,  et  al.  2003,  in  ASP  Conf.  Ser.  294,  Scientific  Frontiers  in  Research 
on  Extrasolar  Planets,  ed.  D.  Deming,  &  S.  Seager  (San  Francisco:  ASP),  427 
Boss,  A.  P.  2006,  ApJ,  643,  501 
Bouchy,  R,  et  al.  2005,  A&A,  444,  L15 

Brown,  T.,  Charbonneau,  D.,  Gilliland,  R.,  Noyes,  R.,  &  Burrows,  A.  2001, 
ApJ,  552,  699 

Bundy,  K.  A.,  &  Marcy,  G.  W.  2000,  PASP,  112,  1421 
Butler,  R.  R,  Marcy,  G.  W.,  Williams,  E.,  McCarthy,  C.,  &  Vogt,  S.  S.  1996, 
PASP,  108,  500 

Butler,  R.  R,  Tinney,  C.  G.,  Marcy,  G.  W.,  Jones,  H.  R.  A.,  Penny,  A.,  &  Apps, 
K.  2001,  ApJ,  555,  410 
Chambers,  J.  E.  2001,  Icarus,  152,  205 

Charbonneau,  D.,  Brown,  T.,  Latham,  D.,  &  Mayor,  M.  2000,  ApJ,  529,  L45 
Charbonneau,  D.,  Brown,  T.  M.,  Burrows,  A.,  &  Laughlin,  G.  2006,  preprint 
(astro-ph/0603376) 

Charbonneau,  D.,  et  al.  2005,  ApJ,  626,  523 


6.  CONCLUSIONS 

Our  study  of  G1  876  has  produced  several  interesting  results. 
We  have  demonstrated  that  transits  of  planets  b  and  c  are  not 
currently  occurring,  and  we  have  thus  modestly  constrained  the 
parameter  space  available  to  observationally  allowed  orbital  con¬ 
figurations  of  the  G1 876  system.  The  absence  of  transits  is  consis¬ 
tent  with  the  dynamical  study  of  Rivera  et  al.  (2005),  who  find  that 
a  coplanar  system  inclination  /  ^  50°  yields  the  best  fit  to  the  ra¬ 
dial  velocity  data.  Our  dynamical  analysis  indicates  that  when  plan¬ 
etary  transits  in  a  strongly  interacting  system  similar  to  G1  876 
are  eventually  observed,  then  a  tremendous  amount  of  informa¬ 
tion  will  be  obtained.  We  look  forward  to  the  discovery  of  the 
first  such  system. 

We  have  also  shown  that  a  globally  distributed  network  of  small- 
telescope  observers  can  effectively  provide  definitive  photometric 
follow-up  of  known  planet-bearing  low-mass  stars.  Our  experi¬ 
ence  with  G1  876  indicates  that  intensive  photometric  monitoring 
campaigns  of  individual,  nearby  M  dwarfs  constitute  a  tractable 
strategy  for  planet  detection.  The  Sun’s  nearest  low-mass  stellar 
neighbors  should  therefore  be  scrutinized  very  carefully  for  transifs 
by  short-period  terrestrial  planets.  A  low-budget  survey  of  this 
type  provides  a  small,  but  nevertheless  nonnegligible  chance  to 
discover  a  traly  habitable  world  ifom  the  ground  and  within  a  year 
or  two.  Such  a  detection  would  be  a  remarkable  and  exciting 
discovery. 
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